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I .   INTRODUCTION 

Inf ragravity  waves  with  typical  periods  of  about  0.5  to 
5  minutes,  often  dominate  wave  runup  on  beaches  during 
storms,  excite  seiches  in  small  harbors,  and  are  believed  to 
be  important  to  nearshore  morphology.  Observed  wave  spectra 
typically  show  strong  variations  in  infragravity  levels 
across  the  continental  shelf  (Figure  1) ,  with  maximum  levels 
close  to  shore  and  weak  infragravity  motions  in  the  open 
ocean.  The  physical  processes  that  control  this  variability 
are  investigated  here  with  data  from  a  cross-shelf  transect 
of  bottom  pressure  recorders . 

Previous  studies  (Okihiro  et  al.,1992;  Elgar  et  al .  , 
1992;  Herbers  et  al . ,  1994,  1995a)  have  shown  that 
infragravity  waves  on  the  continental  shelf  are  a  mix   of 
forced  and  free  waves  (Figure  2a) .  Forced  infragravity  waves 
are  locally  excited  by  non-linear  wave-wave  interactions. 
Two  swell  components  with  slightly  different  frequencies  f 
and  f+Af  excite  a  forced  secondary  wave  with  a  low 
(infragravity)  frequency  Af  (Longuet-Higgins  and   Stewart 
1962;  Hasselmann  1962).  The  observed  energy  levels  of  forced 
infragravity  waves  are  accurately  predicted  by  second-order 
nonlinear  theory  (Herbers  et  al . ,  1994)  but  their 
contribution  to  the  total  infragravity  energy  is  generally 
small  (Herbers  et  al . ,  1995a). 


Free  infragravity  motions  on  the  continental  shelf  are 
either  waves  radiated  from  nearby  beaches,  where  they  are 
generated  by  the  strong  nonlinearities  and/or  wave  breaking 
that  occur  in  very  shallow  water,  or   remotely  generated 
waves  arriving  from  the  open  ocean  (Figure  2a) .  Extensive 
observations  on  both  Pacific  and  Atlantic  continental  shelf 
sites  show  a  strong  correlation  between   infragravity  and 
swell  energy  levels,  suggesting  that  free  infragravity  waves 
are  predominantly  radiated  from  the  beach  [Herbers  et  al . , 
1995a;  and  references  therein] .  Based  on  observations  by 
Munk  (1949)  and  Tucker  (1950),  Longuet-Higgins  and  Stewart 
(19  62)  suggested  that  as  the  incident  swell  are  dissipated 
in  the  surf  zone,  the  associated  forced  infragravity  waves 
are  released  as  free  waves,  reflect  from  the  beach  and 
propagate  seaward  over  the  continental  shelf.  Alternatively, 
Symonds  et  al .  (1982)  suggested  that  slow  oscillations  in 
the  wave  setup,  associated  with  slow  variations  of  the 
breakpoint  location  of  groupy  incident  swell,  can  drive  free 
infragravity  waves.  Although  observations  of  the  directional 
properties  of  infragravity  waves  radiated  from  shore  are  in 
favorable  agreement  with  Longuet-Higgins  and  Stewart's 
hypothesis  (Herbers  et  al  1995b) ,  the  role  of  wave-breaking 
and  associated  set-up  variations  in  the  infragravity  wave 
generation  process  is  still  poorly  understood. 

As  free  infragravity  waves  radiated  from  shore  travel 
into  deeper  water,  they  refract  towards  increasingly  oblique 


propagation  directions  relative  to  the  shoreline  (Figure 
2b) .  Depending  on  the  initial  propagation  direction  in 
shallow  water,  these  waves  may  either  radiate  to  the  open 
ocean  (i.e.,  leaky  waves)  or  reflect  back  towards  the  shore 
from  turning  points  on  the  sloping  beach  and  shelf  (i.e., 
trapped  waves) .  Theoretical  predictions  of  nonlinear 
interactions  in  shallow  water  show  that  directionally  narrow 
swell  incident  on  a  beach  can  drive  a  directionally  broad 
spectrum  of  free  infragravity  waves  that  are  predominantly 
trapped  on  the  shelf  with  relatively  weak  radiation  to  the 
open  ocean  (Herbers  et  al.  1995b). 

Observed  infragravity  energy  levels  on  the  beach,  shelf 
and  in  the  open  ocean  are  qualitatively  consistent  with 
strong  refractive  trapping  and  a  relatively  weak  leaky 
component  (e.g.,  Webb  et  al . ,  1991;  Okihiro  et  al . ,  1992; 
Herbers  et  al . ,  1995a) .  However,  when  local  swell  energy 
levels  are  high,  seaward  propagating  waves  dominate  the 
infragravity  band,  suggesting  that  free  waves  radiated  from 
shore  are   significantly  damped  before  reaching  their 
turning  point  on  the  shelf  (Elgar  et  al . ,  1994;  Herbers  et 
al. ,  1995)  . 

While  infragravity  waves  radiated  from  shore  are 
predominantly  trapped  or  dissipated  on  the  shelf,  a 
relatively  weak  leaky  component  propagates  into  deep  ocean 
basins.  Webb  et  al .  (1991)  show  that  free  infragravity  waves 
observed  in  the  open  ocean  are  the  weak  background  radiation 


from  distant  shores  where  incident  swell  and  infragravity 
wave  energy  levels  are  high.  Distant  (trans-oceanic)  sources 
appear  to  dominate  the  infragravity  wave  field  on  the  shelf 
during  extremely  calm  conditions  when  the  generation  of 
infragravity  waves  on  nearby  beaches  is  relatively  weak 
(Herbers  et  al . ,  1995a, b).  In  contrast  to  the  directional 
broadening  and  trapping  of  waves  radiated  from  shore  (Figure 
2b) ,  remotely  generated  waves  traveling  into  shallow  water 
refract  towards  propagation  directions  that  are 
perpendicular  to  the  shoreline,  causing  a  direct ionally 
narrow,  shoreward  propagating  wave  field  close  to  shore 
(Figure  2c) .  Upon  reflection  from  the  beach  these  waves  may 
propagate  back  to  the  open  ocean  or  (if  the  shelf  is 
irregular)  become  trapped  on  the  shelf.  Observed  directional 
spectra  of  infragravity  waves  on  the  shelf  during  low  wave 
conditions  are  indeed  bimodal  with  narrow  peaks  in  shoreward 
and  seaward  directions,  in  contrast  to  the  broad  directional 
spectra  (with  significant  alongshore  propagation)  typically 
observed  during  high   energy  conditions  (Herbers  et  al . , 
1995a) . 

In  this  study  the  propagation  of  infragravity  waves 
across  the  continental  shelf  is  examined  in  more  detail 
using  measurements  collected  offshore  of  Duck,  North 
Carolina.  Ten  bottom  pressure  recorders  were  deployed  along 
a  cross-shelf  transect  extending  from  the  beach  (6m  depth) 
to  the  shelf  break  (87  m  depth)  for  4  months  during  the  fall 


of  1994.  The  field  data  are  described  in  Chapter  II. 
Estimates  of  the  forced  and  free  wave  contributions  to  the 
infragravity  wave  field  are  presented  in  Chapter  III.  In 
Chapter  IV  the  observed  cross-shelf  variability  of  free 
infragravity  wave  energy  is  examined,  and  the  importance  of 
refractive  trapping  and  damping  effects  is  discussed.  In 
Chapter  V  the  observations  are  compared  to  predictions  of   a 
simple  geometrical  optics  -  based  (WKB)  model  for  the 
propagation  and  trapping  of  long  waves (Herbers  et  al . , 
1995a) .  The  results  are  summarized  in  Chapter  VI.   All 
figures  are  contained  in  the  Appendix. 


II.  FIELD  DATA 

Infragravity  motions  on  the  continental  shelf  are 
investigated  with  field  data  collected  offshore  of  the  U.S. 
Array  Corps  of  Engineers  Field  Research  Facility,  Duck,  North 
Carolina,  as  part  of  the  DUCK94  experiment.  Ten  bottom- 
pressure  sensors  were  deployed  along  a  100-km-long  cross- 
shelf  transect  that  extended  from  the  Duck  beach  to  the 
shelf  break  (Figure  3).  The  three  shallowest  sites  X  (6m 
depth),  A  (12m),  and  B(20m)  span  the  nearly  plane 
(slope=0 . 005)  inner  shelf.  The  next  6  sites  (C-H)  are 
roughly  equally  spaced  on  the  wide,  irregular  (depths  25- 
50m)  outer  shelf,  and  the  deepest  site  I  (depth  87m)  is  on 
the  steep  (slope=0.1)  continental  slope  (Figure  4). 

At  each  site  a  self-contained  battery-powered, 
internal -recording  instrument  package  was  mounted  inside  an 
anchor  with  a  surface  mooring  (sites  A-I)  or  on  a  pipe 
jetted  in  the  beach  (site  X) .  The  instrument  package 
contained  a  Setra  capacitance- type  pressure  transducer,  a 
Tattletale  microprocessor,  and  a  disk  drive  for  data 
storage.  Pressure  data  were  recorded  nearly  continuously 
with  a  2  Hz  sample  rate  during  the  four-month- long  (August- 
November,  1994)  deployment.  Some  malfunctioning  data 
acquisition  systems  were  replaced  with  a  cassette  data  tape 
storage  system  with  a  reduced  sampling  scheme  (a  13 7 -minute- 


long  record  sampled  at  1  Hz  every  3  hours) .  Site  B  suffered 
significant  data  loss  during  the  first  two  months  of  the 
experiment,  and  the  instruments  at  the  shallowest  sites  X 
and  A  failed  during  Hurricane  Gordon  on  November  18 . 
Additionally,  high  resolution  directional  wave  data  were 
collected  in  8m  depth  (1  km  from  shore)  with  a  250m 
aperture,  14-element  array  of  pressure  gauges  (deployed  and 
maintained  by  the  U.S.  Army  Corps  of  Engineers). 

The  analysis  is  based  on  3-hour-long  records  (137 
minutes  at  some  sites  where  longer  records  were  not 
available) .  A  pressure  spectrum  with  0.0005  Hz  resolution 
was  obtained  from  each  data  record  based  on  the  Fourier 
transforms  of  detided  and  tapered  34-minute-long  data 
segments.  The  pressure  spectra  were  subsequently  converted 
to  equivalent  sea  surface  elevation  spectra  using  a  linear 
theory  depth  correction  (the  correction  is  small  at 
infragravity  frequencies) . 

The  strong  variability  of  infragravity  spectral  levels 
across  the  shelf  is  illustrated  in  Figure  5a  with  the 
average  value  of  the  ratio  between  the  spectral  levels 
observed  at  each  site  and  the  spectral  levels  observed  at 
the  deepest  site  I.  While  the  cross-shelf  variations  in 
spectral  levels  is  weak  (less  than  a  factor  of  2)  at  swell 
frequencies  >  0.07  Hz,  a  strong,  consistent  increase  in 
spectral  levels  with  decreasing  water  depth  is  observed  at 
infragravity  frequencies  (0.004-0.04  Hz).  The  infragravity 
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spectral  levels  at  the  shallowest  site  X  are  about  a  factor 
3  0  larger  than  those  observed  at  the  deepest  site  I.  This 
increase  is  qualitatively  consistent  with  both  the 
amplification  of  forced  waves  in  shallow  water  and  the 
refractive  trapping  of  free  waves  radiated  from  shore,  but 
inconsistent  with  the  weak  shoaling  amplification  (°c  h"1/2 
for  a  shelf  with  straight  and  parallel  depth  contours) 
expected  for  remotely  generated  infragravity  waves  arriving 
from  the  open  ocean  (Herbers  et  al. ,1995a;  and  references 
therein) .  While  a  strong  cross-shore  decay  of   infragravity 
energy  levels  is  observed  on  the  inner  shelf  (sites  X-D)  and 
at  the  shelf  break  (sites  H-I) ,  the  energy  variations  are 
more  gradual  over  the  broad  mid- shelf  region  between  sites  D 
and  H  (not  shown) .  These  observations  suggest  that  the 
shelf-wide  variability  of  infragravity  waves  is  primarily 
controlled  by  large  scale  depth  variations  and  not  strongly 
affected  by  local  bottom  irregularities.  The  observed  cross- 
shelf  variations  in  spectral  levels  appear  to  be 
approximately  constant  across  the  infragravity  band.  The 
fluctuations  in  relative  spectral  levels  observed  at  the 
shallowest  sites  X  and  A  are  the  expected  standing  wave 
patterns  associated  with  shoreline  reflections. 

Observed  correlation  coefficients,  between  spectral 
levels  at  infragravity  frequencies  and  the  total  swell 
variance  are  presented  in  Figure  5b.  Swell  variances 
observed  at  the  inner  shelf  site  C  (25m  depth)  are  used 


throughout  this  paper  to  characterize  swell  conditions  on 
the  shelf  because  variations  in  swell  energy  across  the 
shelf  are  generally  small,  the  depth  attenuation  of  swell  in 
2  5m  depth  is   relatively  weak,  and  the  site  C  instrument 
functioned  reliably  throughout  the   experiment.  The  strong 
correlation  observed  at  all  sites  between   infragravity  and 
swell  energy  levels  (correlation  coefficients  squared  are 
above  0.5  in  the  frequency  range  0.004-0.04  Hz)  indicates 
that  the  infragravity  waves  are  predominantly  locally 
generated  (e.g.,  forced  waves  or  waves  radiated  from  nearby 
beaches)  and  that  waves  arriving  from  remote  sources  are 
relatively  unimportant,  consistent  with  earlier  observations 
(Herbers  et  al .  1995a)  .  Note  however,  the  somewhat  weaker 
correlation   observed  at  site  H  and  particularly  site  I  that 
suggests  a  significant  contribution  of  remotely  generated 
waves  near  the  shelf  break. 

At  all  sites  the  correlations  drop  rapidly  at 
frequencies  above  about  0.04  Hz  because  very  low-frequency 
swells  from  distant  storms  are  statistically  independent  of 
the  dominant  (frequencies  >  0.07  Hz)  swells.  The  correlation 
coefficients  are  also  small  below  0.004  Hz,  suggesting  that 
these  low  frequency  motions  are  not  driven  by  swell, 
consistent  with  earlier  observations  (Herbers  et  al . ,  1995a, 
and  references  therein) . 
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III.    FORCED  AND  FREE  WAVE  CONTRIBUTIONS 

The  contributions  of  forced  and  free  waves  to  the 
observed  infragravity  motions  were  estimated  using 
bispectral  analysis.  The  method  is  based  on  the  fact  that 
forced  waves  are  phase-coupled  to  local  swell,  while  free 
waves  do  not  contribute  to  the  bispectrum  (Herbers  et  al . 
1994) .  Analogous  to  the  second-order  energy  density  spectrum 
E(f), 

E(f)df=2E{dP(f)dP(-f)}  (1) 


the  third-order  bispectrum  B(f,Af)  is  defined  as 
(Hasselmann  et  al .  1963)  : 

B(f,Af)dfdAf=2E{dP(f)dP(Af~)dP(-f-Af)  }         (2) 


where  E{  }  is  the  expected  value  and  dP(f)  is  the  Fourier- 
Stieltjes  transform  of  the  pressure  time  series  p(t): 


p(t)  =  fdP(f)exp(27uift)  (3) 
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An  estimate  of  the  contribution  of  forced  waves  to  the  total 
infragravity  variance  (over  the  frequency  range  0.004-0.04 
Hz)  is  obtained  by  integrating  the  bispectrum  B(f,Af)  over 
all  pairs  of  swell  components  (frequencies  f,  f+Af)  with  a 
difference- frequency  (Af)  in  the  infragravity  band  (see 
Herbers  et  al . ,  1994  for  details). 


2fdAffdf  B(f,Af) 
'^  [2J'dAffdf   E(f+Af)E(f)  JdAf  E(Af)]1' 


B«='    '-.  -/'-—  "^     -T-7-      -7-—i  (4) 


The  ratio  between  the  forced  and  total  infragravity 
energies  is  approximately  given  by 


IG,  forced     i  ,-,   i  2  ,  c  » 

E  2jl 


IG,  total 


The  phases  of  the  bispectral  integral  B±i  observed  at 
the  mid-shelf  site  E  (34m  depth)  are  shown  in  Figure  6. 
While  the  imaginary  part  of  BLi    is  small  and  approximately 
randomly  scattered  about  zero,  the  real  parts  include  large 
negative  values.  Although  the  statistical  uncertainty  in 
these  bispectral  estimates  is  considerable,  the  phases  are 
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generally  consistent  with  the   theoretical  value  of  180°  for 
weakly  nonlinear  waves  (Hasselmann,  1962).  Similar  biphases 
were  observed  at  the  other  sites  (not  shown  here)  except  for 
the  shallowest  site  X,  where  small,  but  systematic 
deviations  from  180°  were  observed  that  are  possibly  caused 
by  bottom  slope  effects  (i.e.,  significant  depth  variations 
over  a  wavelength)  or  phase-coupling  between  the  incident 
swell  and  infragravity  waves  reflected  from  the  beach. 

The  estimated  ratios  of  forced  to  total  infragravity 
energy  at  sites  A  (12m  depth) ,  E  (34m)  and  I  (87m)  are  shown 
in  Figure  7 .  The  estimated  forced  fraction  of  the 
infragravity  energy  generally  increases  with  both  increasing 
swell  energy  and  decreasing  water  depth,  similar  to  the 
trends  reported  by  Herbers  et  al .  (1995a).  However,  at  all 
sites  (including  those  not  shown  in  Figure  7)  the  estimated 
forced  wave  contribution  to  the  total  infragravity  variance 
is  usually  less  than  20%  ,  indicating  that  free  waves 
dominate  the  infragravity  band. 
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IV.   FREE  WAVE  VARIABILITY 

Energy  levels  of  free  infragravity  waves  radiated  from 
shore  generally  decrease  with  increasing  distance  from  shore 
owing  to  both  propagation  and  damping  effects.  On  a  gently, 
monotonically  sloping  shelf  with  straight  and  parallel  depth 
contours  h(x),  the  propagation  direction  0  of  a  seaward 
propagating  wave  component  increases  according  to  Snell's 
law,  in  the  shallow  water  approximation: 


sin(0(x))      _   _ 

=constant  (6) 


(gh(x))1/2 


until  the  turning  point  is  reached  where  |0|=7i/2,  and  the 
wave  reflects  back  towards  shore.  Far  from  shore,  the 
combined  effects  of  unshoaling,  refraction  and  trapping, 
result  in  a  directionally  isotropic  wave  field  with  energy- 
levels  that  are  inversely  proportional  to  the  water  depth 
(Herbers  et  al . ,  1995a).  This  h"1  depth  dependence  is 
stronger  than  the  h_1/2  dependence  predicted  for  leaky  waves 
propagating  perpendicular  to  straight  and  parallel  depth 
contours . 

While  the  cross-shelf  energy  variations  owing  to  linear 
propagation  effects  ( (un) shoaling,  refraction  and  trapping) 
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depend  primarily  on  the  directional  properties  of  the  waves 
and  the  shelf  topography,  the  attenuation  owing  to 
dissipation  effects  (e.g.,  bottom  friction)  appears  to  be  a 
strong  function  of  wave  energy  levels  on  the  shelf  (Elgar  et 
al.,  1994;  Herbers  et  al . ,  1995b). 

The  importance  of  both  propagation  and  dissipation 
effects  to  the  cross-shelf  variability  of  free  infragravity 
waves  is  illustrated  in  Figure  8  with  estimates  of  the  free 
wave  attenuation  over  different  parts  of  the  shelf. 
The  ratio  RAD  between  free  wave  energy  estimates  (obtained 
by  subtracting  the  bispectrum-based  forced  wave  energy 
estimates  from  the  total  observed  infragravity  energy,  see 
Chapter  III)  at  inner  shelf  sites  D  (depth  hD=34m)  and  A 
(hA=12m) ,  separated  by  about  30km,  varies  between  0.2  and 
0.6  (Figure  8a)  .  The  majority  of  the  RA  D  values  are  in  range 
0.2-0.5,  roughly  comparable  to  the  theoretical  value 
hA/hD=0.35  for  an  isotropic  trapped  wave  field.  Higher  RA  D 
values  in  the  range  0.5-0.6,  closer  to  the  theoretical  value 
(hA/hD)  1/2=0 . 6  for  leaky  waves,  are   observed  only  on  a  few 
occasions  with  low  swell  energy  levels.  These  observations 
suggest  that  the  infragravity  band  on  the  inner  shelf  is 
dominated  by  waves  that  are  radiated  from  shore  and 
refractively  trapped  by  the  sloping  sea  bed,  and  that 
contributions  of  (leaky)  remotely  generated  waves  are 
significant  only  with  very  low  swell  energy  levels, 
consistent  with  earlier  studies  (Herbers  et  al . ,  1995  a,b) . 
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A  more  quantitative  analysis  of  propagation  and  trapping 
effects  on  the  inner  shelf  is  presented  in  Chapter  V. 

The  ratio  RDF  between  sites  F  and  D,  also  separated  by 
about  3  0  km,  but  in  about  the  same  depth  (hD=34m,  hF=32m) 
varies  between  0.6  and  1.2,  (Figure  8b).  Across  this 
relatively  flat  middle  part  of  the  shelf,  propagation  and 
trapping  effects  are  expected  to  be  weak,  qualitatively 
consistent  with  observed  RDF  values  close  to  1.  However,  a 
systematic  decrease  of  the  observed  RD  F  with  increasing 
swell  variance  suggests  that  waves  propagating  seaward  over 
the  shelf  are  significantly  damped  during  high-energy  swell 
conditions,  consistent  with  results  of  earlier  studies 
(Herbers  et  al . ,  1995b). 

The  ratio  RF H   (Figure  8c)  between  outer  shelf  sites  H 
(hH=49m)  and  F  (hF=32m)  ,  separated  by  about  3  0km,  is 
typically  within  the  range  of  0.5-0.8,  with  again  a 
noticeable  trend  of  decreasing  R  with  increasing  swell 
energy.  The  estimates  suggest  continued  attenuation  of 
infragravity  waves  across  the  outer  shelf  owing   to  both 
propagation  and  damping  processes. 

The  ratio  RH#I  (Figure  8d)  between  energy  levels 
observed  at  sites  I  (hI=87m)  and  H(hH=49m)  situated  near  the 
shelf  break  varies  between  about  0.5-0.7  for  high  energy 
swell  conditions  and  0.6-1.1  for  low  to  moderate  swell 
conditions.  Sites  I  and  H  are  separated  by  only  6  km,  and 
thus,  the  attenuation  owing  to  bottom  friction  is  expected 
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to  be  weak.  The  RHI  estimates,  consistently  larger  than  the 
theoretical  trapped  wave  value  of  0.56,  suggest  there  are 
significant  contributions  of  remotely  generated  (leaky) 
infragravity  waves  at  the  shelf  break.  Although  more 
energetic  waves  radiated  from  shore  dominate  the 
infragravity  wave  field  close  to  shore,  these  waves  are 
attenuated  by  trapping  and  damping,  and  thus,  their  energy 
levels  are  strongly  reduced  at  the  shelf  break.  Relatively 
weak,  remotely  generated  waves  are  usually  not  significant 
on  the  beach  and  shelf  but  may  dominate  the  infragravity 
wave  field  seaward  of  the  shelf  break. 
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V.    PREDICTION  OF  PROPAGATION  AND  TRAPPING  EFFECTS 

To  demonstrate  that  the  observed  strong  decrease  of 
infragravity  energy  across  the  inner  shelf  is  primarily  the 
result  of  refractive  trapping,  the  observations  are  compared 
here  to  a  simple  WKB-geometrical  optics  model  that  describes 
the  transformation  of  a  continuous  spectrum  of  waves  on  a 
gently  sloping  sea  bed  (Herbers  et  al . ,  1995a).  In  the 
shallow  water  WKB  approximation,  the  transformation  of  the 
directional  spectrum  E(6;x)  over  a  monotonically  sloping 
seabed  with  straight  and  parallel  depth  contours  is  given  by 


h 
E(Q;x)  =_7£_E(0s(0f  x)  ;x)  (7) 

h(x)         s  s 


where  0=0  corresponds  to  offshore  propagation,  and  xs  is  a 
location  with  depth  hs  shoreward  of  x  where  the  propagation 
direction  0S  is  given  by  Snell's  law 


h       - 
sin(0  (6;X))  =  (— £-)  2sin0  (8) 

s        h  (x) 


Waves  traveling  seaward  or  shoreward  at  xs  with  propagation 
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directions   0S,    within   the   apertures 


h         1 

|0J  <  arcsin[(-7^-)  2] 
s  h(x) 


(9) 


h         1 
|0  -it |   <  arcsin[(-—^-)  2] 
s  -h(x) 


will  reach  the  offshore  distance  x  while  all  other 
components  are  trapped  between  xs  and  x. 

WKB  predictions  of  the  transformation  of  free 
infragravity  waves  across  the  inner  shelf  were  initialized 
with  estimates  of  the  directional  spectrum  E(6) (integrated 
over  the  infragravity  frequency  range) ,  obtained  from  the  8m 
depth  array  with  a  variational  technique  (Herbers  and  Guza, 
1990).  Two  contrasting  examples   of  directional  spectrum 
estimates  in  8m  depth  are  shown  in  Figure  9 .  The  top  panel 
shows  an  example  of  remotely  generated  waves  with  two  narrow 
peaks  at  about  2  0°  (onshore  propagating  waves  arriving  from 
deep  water)  and  160°  (presumably  the  specular  reflection  of 
these  waves  from  shore) .  The  bottom  panel  shows  a  more 
typical  broad  directional  spectrum  of  waves  radiated  from 
shore  with  significant  energy  at  alongshore  propagation 
directions  (0=90°)  . 

WKB  predictions  of  E(0)  were  made  only  at  sites 
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A(hA=12m)  and  B(hB=20m)  because  this  section  of  the  shelf  is 
monotonically  sloping  with  relatively  small  alongshore  depth 
variations.  The  directional  spectra  predicted  in  12  and  2  0m 
depth  are  similar  in  shape  but  broader  than  the  spectrum 
observed  in  8m  depth  (e.g.,  Figure  9),  and  often  nearly 
isotropic  in  2  0m  depth  (e.g., Figure  9b). 

The  directional  spectra  were  integrated  to  predict  the 
total  free  infragravity  energy  at  sites  A  and  B.  The 
observed  and  predicted  ratio  RA h   between  free  infragravity 
energy  at  sites  B  and  A  are  compared  in  Figure  10  for  the 
period  of  September  22  to  November  17  (when  all  instruments 
were  operational  in  8,  12  and  20m  depth) .  The  observed  and 
predicted  RA B  are  generally  in  good  agreement.  Note  that 
during  days  274-276,  when  local  swell  energy  levels  were 
very  low  and  the  infragravity  band  is  dominated  by  leaky 
(remotely  generated)  waves  (Figure  9a) ,  the  decrease  in 
energy  between  12  and  20m  depth  is  small  (RA  B=0 . 8-0 . 9)  .  This 
calm  period  is  followed  by  a  sudden  increase  in  swell  energy 
and  decrease  in  RAB  (to  about  0.6),  that  indicates  a 
transition  from  remotely  generated  infragravity  waves  to 
waves  radiated  from  shore.  A  similar  transition  occurred 
around  day  314,  but  here  the  predicted  RA B  values  are 
somewhat  larger  than  the  observed  values.  Discrepancies  are 
also  noted  on  days  2  85-287  when  the  observed  attenuation  is 
slightly  stronger  (RArB«0  .  4-0  .  6)  than  the  predicted 
attenuation  (RA  B=0 . 6-0 . 7)  .  These  differences  observed  during 

21 


a  severe  nor'easter  may  be  the  result  of  damping  of  free 
infragravity  waves  between  sites  A  and  B.  However, 
inaccuracies  of  E(0)  estimates  (e.g.,  Herbers  and  Guza, 
1990)  may  also  contribute  significant  errors  in  the  RA  B 
predictions .  Overall  the  comparisons  presented  here  show 
that  the  observed  variations  in  infragravity  energy  across 
the  inner  shelf  are  primarily  the  result  of  linear 
propagation  and  trapping  effects  and  well  described  by  a  WKB 
model . 
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VI.  CONCLUSIONS 

Infragravity  motions  with  nominal  periods  of  0.5-5 
minutes  are  important  to  a  variety  of  nearshore  processes. 
Previous  studies  have  shown  that  infragravity  motions  are  a 
mixture  of  forced  waves  locally  excited  by  non-linear  wave- 
wave  interactions,  free  waves  radiated  from  nearby  beaches, 
and  remotely  generated  free  waves  arriving  from  the  open 
ocean.  The  variability  of  infragravity  waves  on  a  wide 
continental  shelf  is  examined  with  data  from  a  cross-shelf 
transect  of  10  bottom  pressure  recorders,  extending  from  the 
beach  (6m  depth)  at  Duck,  North  Carolina  to  the  shelf  break 
(87m  depth) . 

Bispectral  analysis,  based  on  the  fact  that  forced 
waves  are  phase-coupled  to  local  swell  groups  while  free 
waves  are  un-coupled,  was  used  to  estimate  the  forced  and 
free  wave  contributions  to  the  total  infragravity  energy. 
The  phases  of  the  bispectra  estimates  are  shown  to  be  in 
agreement  with  second  order  non- linear  wave  theory. 

As  in  previous  studies,  the  relative  contribution  of 
forced  waves  to  the  infragravity  energy  increases  with  both 
increasing  swell  energy  and  decreasing  water  depth,  but  is 
usually  less  than  10%  at  all  instrumented  sites.  The 
observed  infragravity  motions  on  the  shelf  are  clearly 
dominated  by  free  waves . 
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The  observed  strong  correlation  between  infragravity 
and  swell  energy  levels,  and  strong  decrease  in  infragravity 
energy  levels  between  the  beach  and  the  shelf  break  indicate 
that  free  waves  radiated  from  shore  generally  dominate  the 
infragravity  wave  field  on  the  shelf. 

The  strong  attenuation  of  infragravity  waves  observed 
across  the  sloping  inner  shelf  is  shown  to  be  caused 
primarily  by  the  refractive  trapping  of  waves  propagating 
into  deeper  water,  and  well  described  by  a  WKB  model.  The 
weaker  attenuation  usually  observed  across  the  relatively 
flat,  but  highly  irregular  outer  shelf  indicates  that  the 
shelf-wide  variability  of  infragravity  waves  is  not  strongly 
affected  by  scattering  from  local  bottom  irregularities. 
However  the  attenuation  of  infragravity  waves  across  the 
wide  outer  shelf  region  increases  with  increasing  swell 
energy.  This  trend  cannot  be  explained  by  linear  propagation 
effects  and  suggests  significant  damping  of  infragravity 
waves  on  the  shelf  during  storms,  consistent  with  earlier 
studies.  Remotely  generated  (leaky)  waves  sometimes  dominate 
infragravity  motions  on  the  shelf  during  periods  of 
extremely  low  swell  energy  levels  when  the  radiation  of 
infragravity  waves  from  nearby  beaches  is  relatively  weak. 
At  the  deepest  instrumented  sites,  weaker   correlations 
between  infragravity  and  swell  energy  levels  and  weaker 
depth  dependence  of  infragravity  energy  levels  are  observed, 
suggesting  that  remotely  generated  waves  are  important 
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seaward  of  the  shelf  break,  possibly  because  the  more 
energetic  waves  radiated  from  shore  are  strongly  attenuated 
across  the  shelf. 
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APPENDIX 

Figure  1.   Example  of  wave  spectra  observed  on  the 
North  Carolina  shelf  (site  locations  are  shown  in  Figure  3). 

Figure  2.   (a)  Sources  of  infragravity  waves,  (b)  Ray- 
trajectories  of  free  infragravity  waves  radiated  from  shore, 
(c)  Ray  trajectories  of  free  infragravity  waves  arriving 
from  a  remote  source.  Dashed  lines  indicate  depth  contours. 

Figure  3.   Bathymetry  of  the  North  Carolina  shelf.  The 
dashed  line  indicates  the  instrumented  transect.  The 
pressure  sensor  sites  are  indicated  by  letters. 

Figure  4.   Cross  section  of  the  instrumented  transect 
(dashed  line  in  Figure  3) .  The  pressure  sensor  sites  are 
indicated  by  letters  (Site  D  is  7.5  km  south  of  the 
transect) . 

Figure  5.   (a)  Average  ratios  between  the  spectral 
levels  observed  at  various  sites  and  the  spectral  levels 
observed  at  the  deepest  site  I  (87m  depth) .  (b)  Observed 
correlation  coefficients  squared  between  the  logarithms  of 
spectral  levels  at  infragravity  frequencies  and  the  total 
swell  variance  at  site  C  (over  the  frequency  range  0.04- 
0.14Hz) . 

Figure  6 .   Imaginary  versus  real  part  of  the  bispectral 
integral  Bi±   observed  at  site  E. 
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Figure  7.   Estimates  of  the  ratio  of  forced 
infragravity  energy  to  the  total  infragravity  energy  at 
sites  A,  E  and  I  versus  swell  energy  (observed  at  site  C) . 

Figure  8.   Observed  ratios  of  the  free  infragravity 
wave  variances  at  sites:  (a)  D  (34m  depth)  and  A  (12m),  (b) 
F  (32m)  and  D  (34m),  (c)  H  (49m)  and  F  (32m)  (d)  I  (87m)  and 
H  (49m) ,  versus  the  swell  variance  (observed  at  site  C) . 
Dashed  lines  indicate  the  theoretical  h"1  energy  variation  of 
an  isotropic  trapped  wave  field  on  a  seabed  with  straight 
and  parallel  depth  contours.  Dashed-dotted  lines  indicate 
the  weaker  theoretical  h"1/2  energy  variation  for  leaky  waves 
propagating  perpendicular  to  straight  and  parallel  depth 
contours . 

Figure  9.  Directional  spectrum  of  infragravity  wave 
energy  E(0)  (arbitrary  units)  observed  in  8m  depth  (solid 
line) ,  predicted  in  12m  depth  (dashed  line)  and  predicted  in 
20m  depth  (dashed-dotted  line)  (a)  2  October  1994  (day  275), 
(b)  4  October  1994  (day  277) . 

Figure  10.   Top  panel:  Observed  swell  variance  (at  site 
C,  units  cm2)  .  Bottom  panel:   Observed  and  predicted  ratios 
RA(B  between  free  infragravity  wave  energy  in  2  0m  and  12m 
depth . 
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